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GOALS VS. ACHIEVEMENTS 



EFFICIENCY <%) 

(sec). 

CATHODE EROSION (Hg/Ct 

GOALS 

50% 

5 000 

10 4 

SELF FIELD MPD 

42% (III Ref. 1) 
30% (A, Ref. 21 

0000 dll' Ref. 1) 
3000 (A, Ref. 2) 

2 X 10 - J <II 2 . Ref. ») 
6 X 10' 4 |Nj, Ref. 7) 
IJ x 10'* (A, Ref. 71 

APPLIED FIELD MPD 

70% (Li, Ref. 31 
70% (II j. Ref. 4) 

50% (Nllj. Ref. 4) 

6800 (Ilj'iur. 5) 

SHOO (l.i, Ref. 3> 
2800 (A, Ref. A) 

3 x 10’* (Hi* 

2 x 1 0* * ( N 1 1 3 * RjL 8) 
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PERFORMANCE FEATURE 

LIMITING EFFECT 

COMMENTS 

TIIRUST EFFICIENCY 

/ 

FROZEN LOSSES 

- IIIOIIEST AT "ONSET” 

• LIMITING 

ionization/kinetic 
ENERGY, (MAY DEPEND" 
GEOMETRY) 


ELECTRODE DROPS 

- FORCE SELF-FIELD MPD 
TO MW POWER 
LESS IMPORTANT WITH 
APPLIED FIELD 

SPECIFIC IMPULSE 

. 

VARIOUS FORMS OF "ONSET" 

1 

. HIGHEST WITH 
LIGHTEST GASES 

LIFE (EROSION) 

ELECTRODE EVAPORATION 

OAS IMPURITIES 
CATHODE MICROARCS 
MASSIVE ARCS AT ONSET 

• THERMAL DESIGN, 
IMPREGNANT 
DISPENSER 

- COMPOSITION 
CONTROL 

- MAY BE IRRELEVANT 
FOR HOT OPERATION 

- ULTIMATE LIMITER 


TIIE MIT PROCRAM 


• SUPPORTED HY AIOSR (JR AMS (19X3 * PRESENT) 

• MAINL Y THEORETICAL WORK. WITH TWO EXCEPTIONS: 

JOINT PROCRAM WITH R & !> ASSOCIATES (IIEIMERDINCER, KILFOYLL) 

JOINT PROGRAM WITH PHILLIPS LAH (CAIDOS) 

• HAS CONCENTRATED ON MODELING FLUID DYNAMICS AND PHYSICS OF SELF-FIELD 

DYNAMICS OF HIGH MAGNETIC REYNOLDS NO. FLOWS EFFECTS OF 
AREA CONTOURING 
EFFECTS OF KINETICS. TRANSPORT 

ANODE DEPLETION AND OTHER HALL EFFECT CONSEQUENCES 
FRICTION. DIFFUSION, HEAT LOSS 
DEVELOPMENT OF MACROSCOPIC INSTABILITIES 

IONIZATION, LOWER IIYHRII) AND ELECTROTHERMAL INDST ABILITIES 
UPPER LEVEL POPULATIONS, INLET EFFECTS 


THRUSTERS: 

I • D MODELS 

1 • D MODELS 

S1AILLUX1 

Kl.UUCS 
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t-D Models', 1 1/2-D Mode If 

\ 

2-0 Models 
(Numerical) 

2-0 Models | Stability | Nadi at ton. 

(Analytical) ; Theory ' kinetics 

i • 

1 

Experimental 

D. Hfimerdinger (Ph.D) 

2BSh33IH 


j 

Contoured 

Channel 

Tic Wins Rmr IMS) 

i 

i 

i 


! 

{ Ionization 
V Stability 


D. Kltfeylc (MS) 

! 


| 

: 

Exit plane 
Spectroscopy 

J.M. Chant; (Ph.D. cant.) 

Ni R h g. ; 

i 

Low 

Interaction 

Asymptotics j 

—■1 


Ell Nlemaed (Ph.D. cand* 

* i 

Physics * 

i 

Nall. 

t-accuratc 



Wmm 


Scan Miller (Ph.D. caad.) 

1 

i 

Transport 

effects 





Jeff Preble (MS) 




Electro- 

thermal 


— 

Eric Sheppard (Ph.D. can< 
Eric Gaidas (Ph.D. cand.) 
M. Marlinet-Sanrhei 

m , ! i lnl ** 

! i | Effect* 


BBBBM 

mm 

1 ! 1 

1 



Onset 

Physics 

II 1 i;li R„ ’ Anode j \ Asymptotics 

Contour- Depletion | 

Electro- 

thermal 

• 


D. Hasting* 

i i 1 

lamer 

Hybrid 




QUASI ONE DIMENSIONAL MODELING 


• BY SACRIFICING GEOMETRICAL DETAIL, EXPLORATION OF A BROAD RANGE OF 
PHYSICAL EFFECTS IS POSSIBLE IN THE CONTEXT OF l-D MODEl-S WITH AREA 
VARIATION 


• SHOWN ARE EXAMPLES OF E. NIEWOOD'S RESULTS ILLUSTRATING 


(a) DEGREE OF AGREEMENT WITH THRUST DATA FROM TWO PRINCETON U. 
THRUSTERS 

(b) RELATIVE IMPORTANCE OF VARIOUS ELECTRON ENERGY SOURCES/SINKS ALONG 
THE LENGTH OF A THRUSTER 























TWO-DIMENSIONAL MODELING - TRANSPORT EFFECTS 


• VISCOUS DRAG IMPORTANT IN SLENDER THRUSTERS 

• VISCOUS DISSIPATION CONTRIIIUTES TO HIGH ION TEMPERATURE 

• DIFFUSION AND HEAT CONDUCTION IMPORTANT AS DAMPING EFFECTS 

RESULTS BELOW FROM S. MILLER'S WORK, FOR D. HEIMERDINGER’S CHANNEL, 
NEGLECTING HALL EFFECT. 

NOTICE BOUNDARY LAYER DEVELOPMENT TO NEAK-FULLY DEVELOPED FLOW. 

LACK OF SYMMETRY IS REAL, AND ARISES FROM ENERGY TRANSPORT BY 
TRANSVERSE CURRENT. 
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Two-Dimensional Viscous MPD Flow 


Magnetic Field 



Max = 0.1 T 
Min = 0.0T 
Inc = 0.004T 


Fluid Velocity 



Max = 5000m/s 
Min = 0 m/s 
Inc = 200m/s 


Gas Temperature 



Max=10000K 
Min = OK 
Inc = 400K 


TWO-DIMENSIONAL MODELING - HALL EFFECT 


• THE HALL EFFECT STRONGLY DISTORTS THE PLASMA FLOWS, AS SHOWN IN THE 2-D 
RESULTS SHOWN NEXT. CONDITIONS ARE 


II = 2 CM. 


L = 10 CM Bo = 0.1 T (I = 30 kA) 
ARGON, m = 4 g/sec 


• NOTICE STRONG AXIAL CURRENT ALONG ANODE. THIS PRODUCES LARGE DISSIPATION 
(SEE T, MAP) AND HIGH IONIZATION FRACTION. PLASMA IS KEPT ELECTROTIIERMALLY 
STABLE BY ELECTRON HEAT CONDUCTION 




VERY STEEP VOLTAGE DROP NEAR ANODE FROM LOCALLY HIGH HALL FIELD. 
POTENTIAL CUT IN NEXT GRAPH. THIS WAS SEEN IN OUR TESTS UNDER SIMIl., 
CONDITIONS (SEE BELOW) 
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[Microscopic Instabilities in MPD Flows 


Microscopic plasma instabilities have been shown to be 
common in many plasma regimes, eg. fusion plasmas, 
ionospheric plasmas. 

In MPD thrusters, current represents a large source of free 
energy, which may drive instabilities. 

Modified Two Stream instability was chosen as a likely 
candidate for importance in MPD. 
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Modified Two 
Stream Instability 


Significant increases in heavy 
species temperature due to 
anomalous heating 

Significant increase in ionization 
fraction due to increased 
dissipation 

Increase in plasma resistivity 
but no macroscopic plasma 
instability 


ORIGINAL PAGE 13 

OF POOR QUALITY 



Conclusions 


• Plasma can evolve to new equilibrium in presence of Mod- 
ified Two Stream Instability, with increased ionization 
fraction and heavy species temperature. 

• Microscopic plasma instabilities cbuld lead to large varia- 
tions in operating voltage and, therefore, efficiency. 

• Plasma instabilities are important in modelling MPD 
flows. 

• Experiments, both existing and, when required, new, 
should be used to ascertain what types of instabilities may 
be excited in MPD flows. 


ELECTROTHERMAL STABILITY THEORY 


UNBOUNDED PLASMA BECOMES STATICALLY UNSTABLE NEAR FULL IONIZATION. 
CONDUCTIVITY - T, 3 '\ SO REGIONS OF HIGHER T, TEND TO CHANNEL CURRENT, 
FURTHER RAISING T, 

______ IC MAQ „cn at PARTIAL IONIZATION BY ENERGY ABSORPTION IN IONIZATION 

PROCESS. SIMILARLY, HEAT DIFFUSION OR ELECTRON-ION PAIR DIFFUSION DAMPEN IT 

FOR SMALL (LESS THAN -2-4 CM) LENGTHS 

WE COUPLED A STANDARD STABILITY ANALYSIS WITH A 1-D MPD MODEL TO PREDICT 
CONDITIONS WHEN 

(a) INSTABILITY WOULD DEVELOP SOMEWHERE (USUALLY AT EXIT) 

(b) GROWTH RATE WOULD EXCEED SOME THRESHOLD 

RESULTS SHOW GOOD AGREEMENT WITH ONSET TRENDS VERSUS 
(a) LF.NOII 

(c) MASS FLOW RATE (THIS DEVIATION FROM I J /ii> SCALING WAS UNEXPLAINED 
BEFORE) 
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SEPARATION OF 'ONSET' AND ANODE DEPLETION 


• FROM OUR COOPERATIVE WORK WITH RiD ASSOCIATES (IIEIMERD1NGER, 1988) 

• USING QUASI - 2D CHANNEL AND 4g/scc. ARGON 

• PROBE AT - 2 mm FROM ANODE DETECTS LARGE AV, AT - 30 KA (CLOSE TO THEORY 
PREDICTION), BUT PLASMA REMAINS "QUIET" 


• AT 60 kA, LARGE, QUASI-PERIODIC VOLTAGE FLUCTUATIONS OCCUR 
VERY CLEAR SEPARATION OF EFFECTS • 


ORICiNAL PAGE S3 
OF POOR QUALITY 
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CURRENT CkA) 


Variation of the Anode Voltage Drop and the Voltage Hash 
aa a Function of the Thmater Currant in the Fully 
Flared Channel for an Argon Maaa Flow Rate 

of 4 g a 


EXIT PLANE SPECTROSCOPIC MEASUREMENTS 



DURING THE SAME TEST SERIES. D. KILFOYI.E USED A 1.26 n». SI'ECTROSCOPE TO 
MEASURE LINE WIDTHS AND LINE INTENSITY RATIOS OF ARGON II AND II LINES III; 
USED AS A DIAGNOSTIC ADDITIVE) 


DATA SHOW HIGH ARGON ION TEMPERATURES (HIGHER THAN T, IN THE ANODE 
REGION). WHICH COULD IMPLY THE PRESENCE OF MICRO-INSTABILITIES 


• DATA ALSO SHOW STRONG ANODE DEPLETION (AT I - 60 kA), IN AGREEMENT >YITII AV A 
DATA 
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Electron Density, crrr^-3 (♦10 1 II *) 


Ti tor FTC ti 9 kV obv|l&| to!U|« 



Rodiol Locotion, 0=cothode 

Radial profile of electron density in fully flared channel for varnmt AAsutti'd 
values of T./T, 


ONSET AS PERFORMANC E LIMITER. PHENOMENA 


AT HIGH CURRENT/LOW MASS FLOW SEVERAL PHENOMENA OCCUR INOT ALWAYS 
SIMULTANEOUSLY WHICH LIMIT *1. Up RANGE 


PHENOMENON 

COMMENT 

! 

(•) 

SHARP RISE IN UNSTEADINESS 

1 

MOST COMMON DEFINITION. 

PLASMA INSTABILITY LOCALIZED 
DOWNSTREAM 

(b) 

i 

INCREASED WALL EROSION 

CLOSELY ASSOCIATED TO U\ 
CURRENT CONCENTRATIONS 

■ 

(C) 

DEVELOPMENT OF LARGE ANODE 
DROP 

NOT ALWAYS PRESENT. 
ALLEVIATED HY ANODE GAS 
INJECTION 

i 

i 

f 

| (d) 

i 

TRANSITION V « I TO V - I 3 

PR OH A III. Y UNRELATED. IIUT 
| HKI.N ASO(TAIEI) W 1 III ONSE 

j 

IIAS 


Al’PROMM \'I E EMPIRICAL CORRELATION: 


AT = 


* if 


I = CURRENT 

m = FLOW RATE 

M = MOLECULAR MASS 

I = accei.etator i. km; hi 

II = IM EN-ELEt l RODE MM AM E 
h = CONS I AM 
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DESIGNATION 


BASIC ASSUMPTIONS 


PREDICTIONS 


COMMENTS 


(a) KQUIPARTITION, OR 
CRITICAL ION1Z. 


VELOCITY 


(b) ANODE DEPLETION 


•ONSET* OCCURS 
WHEN 


eV, = i- 


' \ = 2 0,1 C 


IIALL AXIAL CURRENT - 
FORCES PLASMA AWAY 
FROM ANODE. 

•ONSET WHEN 
<n*>\NODE = W 


Hfrozen = 1/2 

ii _ 2 . /gyr w 

m * p„ V m, i| 


(W « CHANNEL DEPTH) 

(II s INTF.RELKCTRODE 
DISTANCE) 


!L..i _ 945 * k ,T « t\V J 


liH 


16 


o pi 


PROVIDES 

ROUGH 

CORRELATION 

OF 

MOST DAT A 


Al*PRO\|M A I E* 

I.Y SAME 

: EPF.NIlF VTF.S i 

AS DA I A i 

j 

i 

PR ESI Ml* I It )N I 


IS I H A I ANODE 
SHEA I II Nil I 
HR) \R Down 


(C) FULL IONIZATION 


(d) INSTABILITIES 


ONE OF SEVERAL 
ANOMALOUS EVENTS 

OCCURS AS IONIZATION 
ENERGY SINK DISAPPEARS 


SEVERAL PROPOSED: 

. IONIZATION INST All. 

. MICRO INSTABILITIES OK 
1 WO>S I REAM TYPE 

. STATIC ELECTROTHERMAL 


H FROZEN = E IGEOMK I RY) 


li.J- 


fuvy i / her 


W. « DEPTH AT INLET 

A* • TV* II* « THROAT AREA 

; f « Jb- 

Mt,< • VELOCITY FOR 
MOMENTUM m MAGNETIC FORCE 


INSTABILITY THRESHOLD 
IDKPENDS ON TYPE) 

• DOMINANT WAVEIENGIII 

• FREQUENCIES 

• HEATING ETFECIS. FIC. 


• A VARIATION ON **l 


PROVIDES MORE 
DETAIL 


STILL NO 
MECHANISN 


• ELEC 

1 HO| HER 1 1 A I . 
INS MB. 
PROVIDES 
MECHANISM 
FOR H i 

. MICRO. 

INSI AII1I.I1 II S 
MAY EXPLAIN 
HIGH I I 


I 

. Al l. W ROI MILE 
IIY DEI \ll ED 
! probing 


M2 
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[ GEOMETRY EFFECTS ON ONSET | 

BASED ON l-l), VARIABLE AREA MODEL WITH P NEGLECTED. "ONSET" ASSUMED WHEN 

V I - 1- m uj s til c V|/ni| 

TWO CONTOURS: (a) CONSTANT AREA 

(b) CONV. • DIV. (SPACING CHOSEN FOR CONSTANT CURRENT 
DENSITY 

LENGTH MEASURED BY MAGNETIC REYNOLDS N(). BASED ON ALFVEN CRITICAL SPEED: 

(*; * y-^Fl 

TWO MEASURES OF "ONSET” 

(I) NORMALIZED 1 J /m : Y * ^ . IV.f = |i M» A ; 1 

X\ 12 W;,/ ni 

( 2 ) NORMALIZED EXIT VELOCITY: Z = 

'a 

RESULTS SHOW SIGNIFICANT (JAINS IN Hm AND L P BY CONTOURING 

BUT NO GAINS OF P/iii - SHOWING LIMITATIONS OF J* PARAMETER 
GAINS HIGHEST AT LARGE R„ x 


ONSET AT FULL IONIZATION - CONSEQUENCES 


PREBLE’S WORK ON ELECTROTHERMAL INSTABILITY PREDICTS CORRECTLY 
SEVERAL TRENDS. INCLUDING DEVIATIONS FROM l»/rii RULE . 

ELECTROTHERMAL INSTABILITY SEEN TO OCCUR AT a s 0.9 ONLY. 
HOWEVER, ’FULL IONIZATION’ IS NECESSARY FOR INSTABILITY, NOT 
SUFFICIENT. 

(a) GROWTH MAY BE WEAK IN PASSAGE TIME 

(b) IN SMALL CHANNELS OR AT LOW PRESSURES, DIFFUSIVE EFFECTS 
PROVIDE STABILITY 

THEORY STILL TOO CRUDE (LINEAR, CONSTANT BACKGROUND, NO ION 
DYNAMICS...) 

HOWEVER, GIVEN ITS SUCCESS, IT IS INTERESTING TO EXPLORE 
CONSEQUENCES OF ’FULL IONIZATION" MODEL 
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SUMMARY ON SELF-FIELD MPD 


• EFFICIENT ONLY AT HIGH POWER DUE TO LOW VOLTAGE, LARGE ELECTRODE LOSSES. 

• IIIGII POWER OPERATION LIMITED BY "ONSET" ' ' 

• PHYSICS OF ONSET NOT YET CLEAR, BUT IT APPEARS TO DICTATE RATIO OF FROZEN 
LOSS TO KINETIC ENERGY. HOWEVER, THIS RATIO MAY BE CONTROLLED BY DESIGN. 

• ANODE DEPLETION IS SEPARATE LIMITER, ESPECIALLY FOR LONG CHANNELS. 

SHOULD DESIGN FOR COINCIDENT ONSET AND DEPLETION (OR FIND WAY'S TO REDUCE 

AV a \ode) 

• LIFE ISSUES DIFFICULT, BUT PROGRESS IS ENCOURAGING 

• SPECIFIC IMPULSE APPEARS SUFFICIENT IF USING II 2 OR L, 


APPLIED FIELD MPD - THE LOGICAL GROWTH PATH 


! 


i 

i 


• NO TECHNOLOGY' FOR HIGH I*p, COMPACT THRUSTERS IN THE 50 • 100 KW RANGE 


• AF • MPD POORLY UNDERSTOOD, BUT HAS SHOWN POTENTIAL TO FILL THIS ROLE. 
IN ADDITION, NO APPARENT HIGH POWER LIMIT (MAY BECOME SF AT HIGH POWER) 
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(a) INCREASED IMPEDANCE DUE TO U 9 II, VOLTAGE LESSEN IMPACT OF ELECTRODE A V*s - 
SHOULD ALLOW FOR POWER OPERATION. / 


(b) PLASMA ROTATION REDUCES ELECTRODE DAMAGE BY ARCS OR OTHER FAULTS ■ MAY 
ALLOW POST-ONSET OPERATION. 


(c) MAGNETIC CONFINEMENT SHOULD HELP PROTECT WALLS - REDUCE WALL LOSSES. 
LENGTHEN LIFE. 


(d) MAGNETIC NOZZLE SHOULD ALLOW SOME FROZEN LOSS RECOVERY HY EXTERNAL 
EXPANSION. 



(a) ADDED OPERATIONAL COMPLEXITY. BUT SEE RECENT WORK (TAIIARA ET AL., 

ARAKAWA ET AL.) SHOWING POTENTIAL FOR SERIES LOOPS OR PERMANENT MAGNETS. 


(b) INCREASED TESTING DIFFICULTIES (LONG MAGNE TIZED PLUME). BUT LOW POWER 
OPERATION TO COUNTER. 


(c) GREAT PHYSICAL COMPLEXITY . THRUST - PRODUCING MECHANISMS STILL DEBATED. 
SPATIAL-TEMPORAL UNIFORMITY NOT GUARANTEED. REGIMES OF OPERATION 
UNCHARTED. 


1-16 




RECOMMENDATIONS 






REPRODUCE AND VERIFY SELECTED APPLIED FIELD MPD EXPERIMENTS FROM EARLY 
LITERATURE AND/OR FROM AHORAI). 

SUPPORT THEORY/MODELING WORK ON AF THRUSTERS TO EXPLOIT EXISTING 
COMPUTATIONAL CAPABILITIES. 


CONTINUE OUASI-STEADY SF AND AF TESTING TO STUDY DETAILED PLASMA 
MECHANISMS RESPONSIBLE FOR "ONSET" AND OTHER BULK EFFECTS. 


USE 100 - 500 KW STEADY STATE FACILITIES FOR 

(a) STUDIES OF ELECTRODE LIFE AND THERMAL DESIGN FOR BOTH, AF AND SF 
THRUSTERS. 

(b) PERFORMANCE MAPPING AND SYSTEM INTEGRATION FOR AF THRUSTERS. 
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